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Abstract. In this paper, we study the regularity of optimal mappings in Monge's 
mass transfer problem. Using the approximation to Monge's cost function c(x, y) = 
\x — y\ through the costs c £ (x, y) = ^/ e 1 + \x — y| 2 , we consider the optimal mappings 
T e for these costs, and we prove that the eigenvalues of the Jacobian matrix DT e , which 
are all positive, are locally uniformly bounded. By an example we prove that T e is in 
general not uniformly Lipschitz continuous as e — > 0, even if the mass distributions arc 
positive and smooth, and the domains are c-convex. 



1. Introduction 

The Monge mass transfer problem consists in finding an optimal mapping from one 
mass distribution to another one such that the total cost is minimized among all measure 
preserving mappings. This problem was first proposed by Monge [27] and has been 
studied by many authors in the last two hundred years: among the main achievements 
in the 20th century we cite [2T] and |16j . 

In Monge's problem, the cost of moving a mass from point x to point y is proportional 
to the distance \x — y\, namely the cost function is given by 

(1.1) co(x,y) = \x-y\. 

This is a natural cost function. In the last two decades, due to a range of applications, the 
optimal transportation for more general cost functions has been a subject of extensive 
studies. In order to present the framework more precisely, let Q and Q* be two bounded 
domains in the Euclidean space R n , and let / and g be two densities in Q and fl* 
respectively, satisfying the mass balance condition 



(1.2) / f(x)dx = / g(y)dy. 

Let c be a smooth cost function defined on Q x Q*. 
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The problem consists in finding a map T : f2 — > Q* which solves 

min J c(x,T(x))f(x)dx : T # f = g, 

where the last condition reads "the image measure of / through T is g" and means 
J A g{y)dy = f T -i( A \ f(x)dx for all subsets A c Q*. 

The existence and uniqueness of optimal mappings were obtained in [U [TJ [20] if the 
cost function c satisfies 

(A) V (xq, yo) G x fi*, the mappings a; G f2 — >• D y c(x, y Q ) and y G f2* — >■ -D^c^o, y) 
are diffeomorphisms onto their ranges. 

The regularity of optimal mappings was a more complicated issue. For the quadratic 
cost function, it reduces to the regularity of the standard Monge-Ampere equation, 
of which the regularity has been studied by many authors (see for instance |6]). For 
general costs, the regularity was obtained in [26] if the domains satisfy a certain convexity 
condition, /, g are positive and smooth, and the cost function c satisfies the following 
structure condition 

(B) Vx6fl,!/6fi*, and vectors £, rj G lR n with £ _L rj, 

^2 tiiiVM ( c w-s - c p '%- p Cg, rs )c r 'V' / (a;,7/) > /? |£| 2 M 2 , 

i ,j,k,l,p,q,r,s 

where /3q is a positive constant. Loeper [21] showed that the optimal mapping may not 
be continuous if the condition (B) is violated, i.e. when there exist £, rj G M n with £ _L 77 
such that the left hand side is negative. There are many follow-up researches on the 
regularity, in both the Euclidean space [231 E3] an d on manifolds [21 [121 CHI 1221 122] • See 
also [31] for recent development. 

Monge's mass transfer problem is a prototype of the optimal transportation and the 
function (jl.ip is the natural cost function. Therefore the existence and regularity of 
optimal mappings for Monge's problem are of particular interest. However this cost 
function does not satisfy both key conditions, namely the condition (A) for the existence 
and the condition (B) for the a priori estimates. 

The existence of optimal mappings for Monge's problem has been studied by many 
researchers and was finally proved in [HI [32] . Prior to that, the existence was also 
obtained in [16] under some assumptions, and obtained in [30], with a gap fixed in [1]. 
See also [31 [H [10] for the existence of optimal mappings when the norm (11. ip is replaced 
by a more general norm in the Euclidean space. The proofs in [HI [32] are very similar: 
both use the approximation \x — y\ 1+£ — > \x — y\ (e — > 0). The key point is choosing 
an approximation with strictly convex costs of the difference x — y, which satisfy the 
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assumption (A). The optimal mapping for Monge's problem is not unique in general. 
But there is a unique optimal mapping which is monotone on the transfer rays [17] . 

In this paper we study the regularity of optimal mappings in Monge's mass transfer 
problem. As the cost function (II. ip does not satisfy condition (B), the argument in 
[26] does not apply to Monge's problem. Indeed, Monge's problem also admits several 
minimizers T, even if a special one plays an important role: it is the transport map 
which is monotone on each transport ray (see PQ: we will call this map monotone 
optimal trasnport). 

The regularity seems a rather tricky problem and very little is known at the moment. 
Only in the 2 dimensional case, it was proved in [19] that the monotone optimal mapping 
is continuous in the interior of the transfer set (i.e. the union of all transfer rays), under 
the assumptions that the densities f,g are positive, continuous, and with compact, 
convex and disjoint supports. 

Our strategy to attack the regularity in Monge's problem is to establish uniform 
estimates for the optimal mappings with respect to the cost function 



(1.3) c e (x,y) = \/e 2 + \x - y\ 2 

where e G (0, 1] is a constant. The cost function c £ satisfies both conditions (A) and 
(B). Therefore there is a unique optimal mapping T £ associated with c £ , given by 

eDu F 



T £ (x) = x 



1 - Du 



s I 



where u £ is the potential function. By direct computation, u £ satisfies the Monge- Ampere 
equation [26] 

n+2 r 

2 L 
g 



1.4) detwij = ^- 



\Du\ 2 



with 

{»«} = {^ £ 

Under appropriate assumptions, the a priori estimate 

(1.5) sup\D 2 u £ (x)\ <C £ V O' CC n. 

w 

was established in [26], where the upper bound C £ depends on e. Notice that the assump- 
tions involve in particular lower bounds on the densities / and g on their respectives 
domains Cl and Q*. these domains should be c £ — convex w.r.t. each other, which typ- 
ically reduces (if we want to impose it for all e — > 0) to the case of Q C Q*, with Q* 
convex. In particular, this rules out the assumption of [19], since the supports will not 
be disjoint. The case we study is thus completely different form that of [19] . 
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Equation (II Ah is strongly singular as e — > . Note that, due to the small e, a uniform 
bound for D 2 u £ does not mean a uniform estimate for the optimal mapping T £ . Therefore 
we need to work directly on the mapping T £ . 

We wished to prove a uniform bound for DT e , namely the uniform Lipschitz continuity 
of the optimal mapping T £ . By tedious computations, we are able to prove that all the 
eigenvalues of the matrix DT £ , which are all positive, are locally uniformly bounded as 
e —> 0. This is one of the two main results of the paper. Notice that this should bring 
some information on the behavior of DTq, where To is the monotone optimal mapping 
in Monge's problem. Yet, two problems arise: i) the condition on the eigenvalues being 
strongly nonlinear and applied to non-symmetric matrices, it is not easy to pass it to the 
limit, nor to give a meaning to the eigenvalues of DTq (which is a priori a distribution); 
ii) even the fact that the maps T £ do converge to the monotonic optimal transport is not 
that easy if the supports of the measures are not disjoint (which is not the case for us). 

However, as the matrix DT £ is - as we said - not symmetric, the boundedness of the 
eigenvalues of DT £ does not imply the matrix itself is uniformly bounded. Interestingly, 
we find that the matrix DT is not bounded in general. There exist positive and smooth 
/, g such that DTq is unbounded at interior points (here by To we mean the monotonic 
Monge optimal transport, and not the limit of T £ ; however, it is possible to prove (see 
Section 4) that, should DT £ be bounded, then T £ — > T , and hence this implies that DT £ 
cannot be uniformly bounded as e — > 0). This is the second main result of the paper. 

This paper is arranged as follows. In section 2, we state our main estimate, Theorem 
[TJ Section 3 is then devoted to the proof of Theorem [TJ In section 4, we provide 
positive and smooth densities /, g such that the monotonic optimal mapping T is not 
Lipschitz continuous at interior points. We conclude the paper with some remarks and 
perspectives in Section 5. 



Let c = c £ be the cost function given in (II. 3p . The optimal mapping T = T £ : Q — > Q* 
is given by [26] 




2. Uniform a priori estimates 




(2.1) 



T{x) = [D x c{x,-)]- x Du(x), 



where u = u £ is a c-concave potential function. In this and the next sections, we deal 
with the a priori estimates for DT . We will omit the subscript e when no confusions 
arise. 
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From [26] , the potential function u satisfies the fully nonlinear PDE of Monge- Ampere 

type, 

f 

(2.2) det{D 2 x c - D 2 u) = | det D 2 xy c\ in Q. 
For the cost function (11.31) . one has 

(2.3) Dc(x,y) 



\/ e 1 + \x — y\ 2 
Hence by (LP) . 

(2.4) T(x) =x-L(x)Du(x), 
where the function L is given by 

(2.5) L(x) =: 

\J\ — v 

and 

(2.6) v=:\Du\ 2 . 
From (12 .4p and (12. 5p . we can solve 

d 2 (x) 



(2.7) 
and consequently 



e 2 + d 2 (x) 



(2.8) L = ^/e 2 + d 2 (x), 
where 

(2.9) d(x) = \x-T(x)\. 
As in [26] , we denote 

(2.10) A l3 {x) = Dl ixj c(x,T(x)) 

= -^{Sij - DiuDju). 

Then equation (12.21) can be written in the form 

e 2 f 

(2.11) det w { 



lJ L n + 2 goT' 
where 

(2.12) Wij —: — D 2 jU 
is a nonnegative symmetric matrix. 
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We observe from (12. 10D that A^ is positive definite, and the inverse matrix of A^ is 
given by 

(2.13) ,1" L{<y,j I '-jD.uD.u). 
Let us denote 

n 

(2.14) W =: A ^ W *P- 

a,/3=l 

Then we have following uniform estimates: 

Theorem 1. Suppose are bounded domains in M. n (n > 2), f G C 1,l (Q),g G 

C l,1 (VL*), f,g have positive upper and lower bounds, and / ti.jgj) holds. Let u G C 3,1 (f2) be 
a c-concave solution to ( \2.11\) , then we have a priori estimate 

(2.15) W(x) < C, 

where C depends onn, dist(x,dQ), f andg, but is independent of the constants G (0, 1]. 
By (12. 4p (I2.12p and (1 2 . 1 3 [) , it is ready to check that the Jacobian matrix of T is given 

by 

(2.16) Tj = Sij - LjUi - Lu i:j 

L 2 

= - L(uij + —Ui u k u kj ) 

£ k 

k 

Since the matrices {A^} and {wij} are positive, then DT is diagonalizable, and its 
eigenvalues Ai ) • • • , A n of Jacobian DT are positive, and Yli=i ^« = W . So if W is 
bounded, one immediately sees that all the eigenvalues of DT are bounded from above 
and below. We therefore have 

Corollary 1. Under the assumptions of TheoremUi we have for any Q' CC Q, 

(2.17) C' 1 < minAi < maxA* < C in SI', 

i i 

where C depends on n, dist(Q! , <9f2) , / and g, but is independent of e G (0, 1]. 
In view of f l2TT3|) and (12H| . one finds that 



W = w a + — u i u i w ir 



hi 
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Hence we obtain 

L 22 W H - C 

By (12.101) and (I2.12p . we also obtain the estimate for D 2 u. 

Corollary 2. Under the assumptions of TheoremUi we have for any fi' CC O, 

(2.18) \D 2 u\ < C/L in ft', 

where C depends on n, dist(Q' , <9ft) , f and g, but is independent of e G (0, 1]. 

Corollary 2 morally gives a C 1 ' 1 estimate for the potential function uq = lim £ ^o u on 
the set 

E s = {x ett: \T{x) - x\ > 5 > 0} . 

This recovers a well-known result which reads "the potential is C 1 ' 1 in the interior of 
transport rays", which was also used by [1] in order to prove the countable Lipschitz 
property of the direction of Du. At a point xq with v(xq) > 0, denote 

Du (x ) 

\Du {xq)Y 

and let £ Q be unit vectors such that {u, C, a } a=1 n _ 1 are orthonormal. We denote 



By (12TT6D and (l2~T3l) . 

(2.19) D u (u, T) = T^ = J2^ lk ^ 3 



L Y1 ( 1 + V kW kj Vj 



i,j,k 

L 3 



Similarly, 

(2.20) (r, t) = t|; = lJ2 

j,k 

Noticing that {wij} is positive definite, it is clear from (I2.19P and (I2.20p that T" and 



TL are positive. Recall that 



W = T» + Y J T i 



a=l 
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Hence by ( 12. 15ft we obtain 

Corollary 3. Under the assumptions of TheoremUi we have for any fi' CC fi, 

r t v v < c 

(2.21) <C m ' 

where C depends on n, dist(fl',dfl), f and g, but is independent of e G (0, 1]. 

At the limit, this corresponds to saying (even if what we state here is not rigorous) 
that the limit mapping T is Lipschitz continuous in the direction of transfer rays, and 
for any unit vector £ perpendicular to the transfer rays, (£, T ) is Lipschitz continuous 
in the ^-direction. The Lipschitz continuity along transport rays is not surprising, since 
we are doing one-dimensional optimal transport between two measures with upper and 
lower bounds on their densities; yet, the densities of the one-dimensional problem along 
each ray are affected by a Jacobian factor (due to the decomposition of / and g along 
rays), and this makes this Lipschitz result not completely evident. In section 4, we will 
construct an example to show that the component (uq,Tq) is in general not Lipschitz 
continuous in £, even though the mass distributions are positive and smooth, where v 
is a direction of transfer rays and £-Lz/ . 

In Theorem [H we assume that u G C 3 ' 1 . This assumption is not needed if Q C Q* and 
Q* is a bounded convex domain in M. n , as it implies that Q* is c*-convex with respect to 
Q and by approximation, and the condition u G C 3,1 is always satisfied, see [26] . 



3. Proof of Theorem 1 

To prove Theorem (TJ we introduce the auxiliary function 

(3.1) H(x) = r]{x)W{x), 

where 7] is a cut-off function. Suppose that H attains its maximum at some point xq. To 
prove that H(xo) is uniformly bounded in e, the computation is rather complicated. We 
find the computation can be made a little simpler if we first make a linear transformation 
such that 

(3.2) Aij(xo) = 5ij, 
and then make a rotation of coordinates such that 

(3.3) Wij(x ) = diag {X u A n } . 
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It is well-known that Ay, Wij are tensors [22]. An advantage of working on tensors is 
that one may choose a particular coordinate system to simplify the computation. As 
we only made a linear transform on the Euclidean space R n , the Riemannian curvature 
tensor under the metric er^ vanishes, which allows us to exchange the derivatives freely. 
In the following we will use D to denote the normal derivatives in W 1 and V to denote 
covariant derivatives under the metric a. 

Suppose the linear transformation is given by y — P~ x x (i.e. = Yl PikUk) such that 
P T AP = I is the unit matrix at x . Then by ( I2TTU1) and (12121) . 

Aij = J2 A ki P kiPij = ( pTAP ) ir 

k,l 

= wkiPkjPij = (p t wp).. , 

k,l 

where bar denotes quantities in the ^-coordinates. 

Denote {a.^} = P T P, and {a^} = (P T P)~ 1 . Then by fl2TT0D and (12151) . 

(3.4) <),, = A ij ^(P T AP). j 
and 

(3.5) = A« = {P- l A-\P T Y l ) i3 

= L{o'i ■ '—uu'). 

where u l =: Ylk^^k- Note that by (12. 7p and (12. 8p . v and L are invariant under the 
coordinate transformation, and 

(3.6) v = a^UjUj = '^^u l u i < 1, 

(3.7) e 2 < L 2 < C. 

For simplicity we will omit the bar below. In view of (13. 4p . we have, at Xq, 

(3.8) o"jj = LSy + ViuVjU. 
By (13.51) . we have, at Xo, 
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where Ui = Vjix. By (12.51) . it follows that 



Tn 1 T, 3 

(3.9) Ul = ^- = V 

1 — v e 2 



2 

Hence u l = j^Ut and by (13. 5p 



(3.10) a« = I^-^u,^). 

Formulas (13.81) . (13.91) and (13. 10|) will be repeatedly used in our calculation below. With- 
out loss of generality, we may also assume 

(3.11) Ax>A 2 >--->A n 
at xq. 

Since Xq is the maximum point, we have 

77- W- 

(3.12) = V l \ogH(x ) = ± + ^ 



T) W W 

> V ?.log^(x ) = ^-25f + 

■q r] z W 



as a matrix, where subscripts i,j on the R.H.S. denote covariant derivatives in the metric 
a. We thus obtain, at x , 

(3.13) > W/^u^(^-2^) + ^7i/W 4J , 



where w tJ is the inverse of Wij. 



Differentiating ( 12. lip gives 
(3-14) 

(3.15) ^2w V W ij]ab = } y W^W^Wjj.qWst-b + (fab, 

where tp is given by 

/ e 2 f 

(3.16) ip = log 



,L n+2 goT, 

In our computation we use the notation Wy-k = VkWij, Wij-ki = V/VfcWjj, A^ = VfcAj 
and Aij-ki = ViV k Aij. 

To estimate the term ^w^Wij in (13.131) . we first prove the following lemma. 
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Lemma 1. We have 

L 2 1 

Aj;k = -^AijU h W hk + —{UjW ik + UiWjk) 
6 Li 

——(AijUk + AikUj + UiAjk), 

L 2 1 

Aa-p - Aip-i = —(Aurfwtp - Aipi/wu) + y(wipUi - waup). 

E ±j 

(we use the summation convention u h Whk — 2^ih uhw hk)- 

Proof. Recall that v = a^UiUj = ^2 u l Ui. Therefore 

dL 1L 3 

dv ~ 21 2 ' 
v k = 2u h u hk . 

By (13 .4p . Ay = — UiUj). Differentiating, we get 

(3.17) Aij-k = ~j2~^ Vk ( a V ~ UiU i) + -j^(- U ikUj - UiU jk ) 

L? 1 
= — -AijU h u hk + —(-UikUj - UiU jk ) 

E Li 

L 2 1 

= —AijU h w hk + —{wikUj + UiWjk) 

£ Li 

— — (AijUk + AikUj + UiAjk). 
The second formula follows from (I3.17P immediately. 

Differentiating (12.141) twice and using A^ (x ) = 

(3.18) Y, "• i;lr " = «•""•«,:,, - ■>Y," JJ -\« l:i>r n | :j - ^ W ij A 

+2 ^ w l] A a k-iA pk -jW a p 

We have by flBTTBD 

(3.19) V. . w ij w aa .ij = V w ij A aa .ij - V 

^ y w ^u)jj- aa -\- ^ w ^ (A aa> ij Aij- aa ^j 



ijVJ aa . 



w lJ u 



aaij 
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By the first formula in Lemma [1] 

L 2 1 

E. . W lj 'AafriWafrj = — } W aa -iU l - — V] W tJ UjW aa]i 

2 2 v 

+ — UpW a p. a - — 2^ W^UpWjp.i. 

By d33D, it follows 

L 2 1 

z — 'i,3,oi,P £ z — ' -L z — ' 

2 2 

+ — ^ ^(A^q - v4 aa;/ g) - — ^ 

2 

By (13.141) and the second formula in Lemma [U, we then obtain 

L 2 1 



where 



Recalling fl3D and f l3~6|) . 

(3.20) < u'ui < y^u'ui < 1 
for any given z. Hence 

L 2 1 

(3.21) w 11 ' Aap.jWap-i < 3^y^w aa -iU % -—y^w ll UiW 

Here and below we use Q to denote quantities satisfying 

Q <c(l + — + W 2 + -WW) . 

\ n n / 
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Inserting (133311 and (pOIj) into (13351) . we obtain 

(3.22) ^W^Wij > ^^(Awitf- A«;a a ) -J^^-AaajiiWaa 



+ 2^^ + — z^^ M 

To proceed further, we need the following lemma. 

Lemma 2. We have 

} w ii ( y A aa . ii - A ii]aa ) > —W w^UiW^i 

L 2 L 2 
-^W ^ w^;^ + (n + 2) — 2^ ¥>X 

and 

Proof. In view of Lemma (U 

(3.23) A aa yi = —A aa u Ufa —u a u. 



e 2 ™ m L 



at ■ 



By differentiating ( I3.23p . 

L 2 h L 2 

(3.24) A aa -a — A aa U Uhii 7rA aa 0~ Uti^hi 

L 2 / L A t h 

n'A aa .iU 2 —A aa u Unit Ufa 
e 2 e 4 

2 2 2 2L l 

L L e z 

Plugging (13T23|) into (13331) . we infer that 

r2 r2 
4 A h a ht 

L4 ■ ..fc 2 



HQ: 



By (JHUDD, 



^ uu — — t,t — — ill y 

S Lj 

2 2 h 
L e z 



L? L L 4 



14 

Hence 



QI-RUI LI, FILIPPO SANTAMBROGIO, AND XU-JIA WANG 



A 



^ 2 Ah L a 2 

6 £ 

2 2 2 4L • 

L L e z 

j2 t2 t 



2 2 



2 2 4L 
L e z 



Employing (13.231) again, it follows 

L 2 



(3.25) 



.4 



Hence 

(3.26) V w ll A aa .i 

< * i rv 



^4 2 _ 2 2 4L 

e z L e z 

H 7 A aOL AnUhtu u + 2— A aa uaUiU 

L a 4 

e 2 L 2 



L 2 L 2 

—W fhU h + 2— ^ iu«uVi 



L 

'72 



H^W uaUiU 1 + 2W WaUaUiU 



Recalling (13~^0|1 . 



and 



}<Q, 



WW Y u a u i ui + 2W WaUaUiU 1 < - 1 WW + 2 W u^n 1 J> + Q, 



the second inequality of Lemma [2] follows from (13.261) . 
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From (EZSD it follows that 



A A ^ A h ^ 

£ Li 

£ 2 L 

-L. 2 4L • 

e z e z 
L 2 4L 



Hence 



(3.27) ^ ^ U> (^4aa;M — ^4ii;a«) — ^ ^ ] 



L 2 2 

— W 2^ Waa;hU h ~ — ^ W n U i W OLa - i 



e 2 — 



Since 



and 



the first inequality of Lemma [2] follows from f!3.27j) . □ 
In view of Lemma [21 (I3.22p can be rewritten in the form 

(3.28) J2 wijWi i ^ ^Wj2< + ^(Wj2 w ii u i ui + 2 J2 w ii UiU ^ 

-^(W + 6) V* + ^ - ^<2, 

where 

2L 2 L 2 

(3.29) 3f^=: -— ^u;tfuVi + (w + 6 -^)^E^ + X ^ 



By (gUP , we have 



(3.30) X = X ^"i^" ~ 
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It follows from ( EZSj) that 

rVk 



(3.31) 



a 



L 2 2 
-W— -Wu kk u k - —w kk u kk u k 

Tj 6 L 

-W 7 ^ + %Ww kk u k - ^-Wu k 



21/ 



wLu k 



2L? 



-W kk U 



Hence, by (E3D and (Ig^O]) 
(3.32) 



Therefore, by inserting (I3.32p into (13. 28 p . we find that (I3.13P can be written as 

L 



(3.33) 



?7 e^ 



Without loss of generality, we may assume the cut-off function r] satisfies \Dr]\ 2 < 
(otherwise we may replace rj by r) 2 ) and \D 2 r]\ < C . Hence it follows 



w 



13 



■Vij 



2^) + ^(W + 6)^^ 



> -C 



V 

> -±Q, 



\D 2 r)\ 



\D V [ 



L 2 W 



(W+l) 



where (13. 8 j) is used in the last inequality. Therefore (I3.33P can be written as 
(3.34) >^W^^1 ^ 

Lemma 3. M^e /iave, at x , 



L 



Q. 



Proof. Recalling (12. 6p . we have 
(3.35) 



2u h u ha . 
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Differentiating ( 13. 16ft gives 

(3.36) = (/I + 2) — li Uha H 7" ~ — , 

(3.37) y a p = -2{n + 2)—u h u t u ha u t p- y -—^—v a p+ J -^f 

£ A E J 

Up V 2 a p{9 °T) + V a (go T) Vpig o T) 



P 9°T goT goT 

Inserting (I3.36j) and (I3.37|) into (13.291 ). we obtain 

(3.38) ^ > 2— 2^ w «— -(n + Q-W)— ^ -^ T - 

v^ VygoT) (n + 2) L 2 ^ L 

Differentiating (I3.35p . we obtain 

^2 Vaa = 2 ^ ^UtaUha + 2 "^oa* 

= 2 2^ cr th u ta u ha + 2 2^ u h A aa;h - 2 ^ u h w aa . h . 
By (ED, (T3~T0D and (T3~20l) . 

^UtaUha = J ~ U^Ui) < -Q. 

From ( )3~23l) . 

Ej n + 2 4 1 

w A aa] i = — — Uj < ~j^Q- 

Also, by (I3T3T]) . 

-J>^W = wY, — -i(W-2)J2vk k u k u k 
+jWv - -Jjtujjtuu* 

Therefore we have 

(3.39) J2v aa <-Q. 
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It then follows from ( 13. 38ft 

2L 2 ^ V t {goT) 



(3.40) ^ > -rE 



-u 



e 2 ^—^ goT 



—in 



e 9° ± 

VL(soT) L 



' JO J e z 

Now we compute V Q (g o T) and oT). By ( 12.161) we have 

(3.41) V a ((? o T) = 3 fe T Q fe = g a w aa . 

By differentiating (12.161) . we have 

= ^ 9kA kl wi a]a + ^ 9 aa w 2 aa - ^ 9kA ka;a w aa . 
Recalling that A = 5m at xo, we have 

A Wi a - a Ulkot;a ^aa;k A ka - a ^aajfc' 

By (J330D, 

^ V^ Q (# O T) = ^ 9kW aa - k + ^ ^(Afcaja - A aa . k ) 
+ ^ 9aa.W 2 aa - ^2 9kA ka ■ a W aa 

= ~w y^fl%— + g aa wi a + ^ 9k{A ka . a - A aa . k ) 

^ ^ gk(A aa - k ^4/ca;a)^aa- 

Using the second formula in Lemma [TJ we get 



(3.42) ^VL(^oT) = + + 
Inserting (l3Hj) and (I3T42]) into (IBTiUj) . we then obtain 

„ o r 2 

(3.43) (ffo^ > wY J 9k V - + - r J2<9 k u k 

Tj £ 

^2" 2^ WkkU 9k ~ 1^ 9aaW aa 
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By (E3}, 

(3.44) J>3* = ^B^ + F^*! 

'/ fc 

We have 

Y,° ii 9irii = iPg,Dri)<C, 
where D is the normal derivative in W 1 and (•, •) denotes the standard Euclidean metric. 

Similarly, Yli a ^ u i9j = 12 u ^9ji Yli (T ^ u i r lj = ^2 u ^Vj an d 12 a ^9i9j are a ^ bounded by a 
universal constant C. Hence from (I3.44p . 

3.45 >9k—> — o — ■ 

Employing ([33]) and ( 13T201) . 

(3.46) (u fc ) 2 = < ^, 

for any given k. Using (13.101) then (13. 9p . we have 

It implies 

(3-47) = *B*<Wi ■ ^B"'-'/<) 2 

T 3 

< c L - 

e z 

By fl3~46l) and (l3~47j) it follows that |# fc u fc | < C. Hence 

(3.48) E^ fc ^>-^ 2 
and 

(3.49) J2w kk u k g k <CW. 
Moreover, in view of (13. 8p and (I3.9p . 

L 3 L 

au = L + -irUiU 1 < C^r, 

e 2 e z 

for any given i. Consequently, 

(3.50) Baa^L " I^B"^ < C^ 2 ' 
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By virtue of (B35D , (E3HD, ((3321) and (133U1) . we obtain from (EH that 

L 

> -4c 

e 2 

This completes the proof. □ 
By Lemma [3] and (I3.34J) . we get, at x , 

L W 2 L f W o W 
> C-(l + — + W 2 + — 

e 2 n e 1 \ r/ r/ 

Multiplying nifL to both sides of the above inequality, we obtain 

(3.51) > ^rW(H 2 -CH)-C^r(l + H 2 ) 

e z e z 



C^rWH 2 -C^r(l + H 2 ). 



e 2 e 2 



Note that by (ETTB . 

(3.52) ^>Y.i>-(Uif>-^>c( W 



k>2 K k>2 



V n 



n-l 



where C is independent of e. Hence from ( I3.5ip we get 

(3.53) >^H 2+ ^-C^(l + H 2 ). 

Hence H < C at xq and this completes the proof of Theorem [TJ 



4. A COUNTEREXAMPLE TO THE LlPSCHITZ REGULARITY 

In the last section we proved that the eigenvalues of DT £ are uniformly bounded. In 
this section we give an example to show that the T E is not uniformly Lipschitz continuous 
for small e > 0, i.e., the matrix DT e is not uniformly bounded, even though the densities 
/ and g are smooth and positive, and the domain Q* is c-convex with respect to Q. 

Our counterexample will be obtained by finding a choice of / and g such that the 
monotonic optimal transport T between them is not Lipschitz continuous. Even if we 
said that the convergence T e — > T is not straightforward, we can prove that a uniform 
Lipschitz bound on T £ would imply such a convergence, and hence the same bound on 
Tq. Hence, if Tq is not Lipschitz, then T e cannot be uniformly Lipschitz. 
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Lemma 4. Suppose that the sequence of transports T £ is uniformly Lipschitz. Then the 
whole family T £ converges uniformly as e — > to the unique monotonic optimal transport 
for the cost \x — y\, which will be Lipschitz with the same Lipschitz constant. 

Proof. By Ascoli-Arzela's Theorem, the uniform Lipschitz bound implies the existence 
of a uniform limit up to subsequences. Obviously this limit map T will be optimal for 
the limit problem, i.e. the Monge problem for cost c(x,y) = \x — y\ and will share the 
same Lipschitz constant as T £ . 

We only need to prove that T is monotonic along transport rays. Take L £ (x) = 
\J E 2 + \T e {x) — x\ 2 : these maps are also uniformly Lipschitz and converge uniformly to 
L(x) = \T e (x) — x\. Let us denote by u £ the potentials for the approximated problems 
and by u the potential for the limit problem. Due to the uniqueness of the Kantorovich 
potential u, since all the functions u £ are 1— Lipschitz, we have u £ — > u uniformly. 
Moreover, Du £ — ^ Du and the convergence is actually strong (in L 2 , for instance) if 
restricted to the set Tu = {\Du\ = 1} (as a consequence of \Du £ \ < 1, which implies 
that we also have J Tu \Du £ \ 2 —> f T \Du\ 2 : this turns weak convergence into strong, and 
hence also implies pointwise, convergence). 

The monotonicity of T is proven if one proves DL ■ Du < 1, since the direction of the 
transport rays is that of —Du. This inequality is needed on the set of interior points of 
transport rays, which are exactly points where \Du\ = 1. On these points we can use the 
weak convergence DL £ — ^ DL (weakly-* in L°°) and the strong convergence Du £ — > Du, 
which means that it is enough to get DL £ ■ Du £ < 1, and then pass the inequality to the 
limit. This is the point where we use the uniform Lipschitz bound on T £ : without such 
a bound we could not have the suitable weak convergence of DL £ . 

In order to estimate DL £ , we use (12.81) and (12 .4p . We come back to the notation 
without the index e, and write DL, thus getting 

DL ■ Du = -D,u (Tj - <%) DjU = LD,u D 2 -u D 3 u + \Du\ 2 DL ■ Du. 

Then, we use (12. 10)) and (12 . 12)) and the positivity of the matrix Wij, to get 

LDiU DfjU DjU < \Du\ 2 {\ - \Du\ 2 ). 

This implies 

(1 - \Du\ 2 ) DL-Du< \Du\ 2 (1 - |£>w| 2 ), 

which provides DL ■ Du < \Du\ 2 < 1 (notice that, for fixed e > 0, the norm of the 
gradient \Du\ is strictly less than 1, which allows to divide by 1 — \Du\ 2 ). □ 
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To construct the counterexample where To is not Lipschitz, our idea is as follows. Let 

(4.1) 4 = {(x,y) in R 2 | y = \fa (x + 2 + a) ,x £ [—2 — a, 1] } 

be a family of line segments £ a , where a G [0, 1]. It is clear that the segments l a do not 
intersect with each other and U ae [ 0i i]^ a = Aabc, where Aabc denotes the triangle with 
vertices A = (-3,0), 5 = (1,4) and C = (1,0). Let 

/ = 1, 

1 

9 = l + jX + r)(y) 

be two densities on Aabc- We first show that there exists a smooth positive function rj 
such that /, g satisfy the mass balance condition 

(4.2) f f=f g, for all a e [0,1]. 

J ^PaCQa J A PaCQa 

Here P a = (—2 — a, 0) and Q a = (1, (3 + a) y/a) are the endpoints of £ a . We then prove 
that there is a Lipschitz function u, which is the potential function to Monge's problem 
in Aabc, with the densities f,g given above. By (14. 2[) we can construct a measure 
preserving mapping To, which pushes the density / to the density g, with {£ a } as its 
transfer rays. Using the potential u and the duality we show that To is the optimal 
mapping of Monge's problem. By reflection in the a;-axis, we extend T to the triangle 
A abb 1 i where B' = (1, —4) is the reflection of B. Then T is not Lipschitz at the interior 
point (-2,0). 

Lemma 5. There exists a smooth positive function rj, such that holds. This 

function satisfies rj(y) = 0(y 2 ) as y — > 0. 

Proof. By direct computation, 

/ / = 7^(3 + a) 2 , 

J^PaCQa z 

/ g = / / i + + v (y)) dy dx 

J&P„CQ n J -2-a JO V 4 / 



,1 

(3 + a) 2 (12 -a) + / / r]{y)dydx. 



24 J -2-a JO 

In order that ( 14. 2 p holds, we need 

j rl /•y/a(x+2+a) 

(4.3) a 3 / 2 (3 + a) 2 = / / V (y)dydx. 

^ 4 J -2-a JO 
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Differentiating ( 14. 3 j) with respect to a, we have 



a 
24 



(9 + 7a)(3 + a)= / (x + 2 + 3a) r] {yfa (x + 2 + a)) dx 

J-2-a 



which is equivalent to 

(4.4) — (9 + 7a) (3 + a) = / (t + 2 av ^) 77 (t) dt. 

In order to find 77 satisfying (14. 3 P for all a G [0, 1], we only need to solve (14. 4p . since the 
equality in (14. 3p is true for a = 0. 

Let us introduce 



(4.5) y = (3 + a) y/a. 

It is clear that y is a strictly increasing function of a. Let a(y) = 0(y 2 ) be the inverse 
function of ( 14. 5ft . Differentiating ( 14. 4ft in y and using a y = , we obtain 

^ (27 + 45a + 14a 2 ) = ^7^- V GO + f »? (*) ^ 
36 2i/a J 

Taking derivative again, we obtain 

(4.6) ^( l/ ) + i^M^( y ) = Mp ( ( y )) > 

y 

where 

(5a — 1) (3 + a) 

!(°) = — Qn , ^2 — > 
3 (1 + a) 

27 + 135a + 70a 2 



162 (1 + ay (3 + a) 
Solving ( I4.6p . one finds an explicit formula for rj: 

(4.7) r](y) = j\p (a (t)) exp f- jf * q ^ a ^ drj dt. 



It is clear that 



Hence 



-1 <q{a{y)) < if \y\ « 1. 
0<v(y) < C jf tex p(^ ^rj a 7 /: 
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From (14.71) it follows that 



"y ( i r a (v) 5 a — l 
(4.8) Tj(y) = I tp(a(t))exp -- / — -da \ dt 



\ 2 Ja{t) a (1 + a) 
y/a /' a 3(s + l)(27+ 135s + 70s 2 ) 



324(a + l) 3 i v 7 ^ 
a (10a 3 + 41a 2 + 54a + 27) 
54(a + l) 3 

In the last two equalities, a is the function of y determined by (14.51) . Therefore r\ is 
positive and smooth and satisfies the required conditions. □ 

Remark 1. From (f^.5| ) ; we can explicitly write 

1 

= h(y) 



h{y) 

where 



h{y) = \j\l\y A + y 2 + \y 2 + i- 

It is clear that a (y) is a smooth even function. 

Lemma 6. There exists a function u : A^bc ~^ K satisfying 

\u (p) - u (q) | < \p - q\, V p,q G A ABC , 
and equality holds if and only if both p and q lie on a common segment l a . 

Proof. We will construct a function u : Aabc ~^ which decreases linearly along all 

4- 

For (x, y) G Aabc, let a = a (x, y) be the solution of the equation 

(4.9) y = y/a(a + 2 + x) . 

Hence (x, y) G £ a . Differentiating ( 14. 9 p with respect to x and y respectively, we get 

(4.10) = ^L y + y /E(a x + l) 
and 

(4.11) l = %j + Va~a y . 
It follows by (T4TT01) and flUD]) that 

(4.12) aj, + = 

provided a (x, y) ^ 0. 
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On the other hand, for (x, y) G Aabc, the direction vector of the line segment £ a 
passing through (x, y) is given by 

1, y/a (x,y) 



(4.13) v(x,y) = (vi,u 2 ) = 

V 1 + a \ x ^y) 

Hence, by 

(4.14) - a^ 2 = — — — „ 

2(l + a) 3/2 V v 7 ", 

provided a (x, y) ^ 0. 

Fix a point P = (-2, 1). Let 

7(*) = 7x(t) 

= (i(x + 2)-2,l-t(l-y)), 
t G [0, 1]. Then 7 is the segment joining P and X = (2, y) G Aabc- Set 

u (x, y) = (x + 2) [ v x (7 it)) dt + {y-l) [ v 2 (7 (t)) dt. 
Jo «/o 

We claim that u satisfies 

(4.15) Du (x, y) = v (x, y) on all segments £ a . 

Indeed, for any point X = (x , y ) G Aabc with a (x , y ) 7^ 0, by (14.140 we have 

(4.16) u x (x ,y ) = [ v 1 (y (t))dt+ (x + 2) f td x v x { lo {t))dt 



+ (2/o-l) [ td x v 2 ( lo (t))dt 
Jo 

vi(jo(t))dt+ / t— 1/1 (70 (*)) dt 
1 d 

37 Oi (7o (<))) d< = fi (^0, 2/o) , 
dt 

where 70 = 7x an d we used d x u 2 = d y v\. Similarly, we have 

u y (x ,yo) = "2 (xo,yo) ■ 
Taking limit, we see that (14.150 also holds on the segment £ a \ a=0 . 
As v is a unit vector, hence from (14.151) we have 

\u (p)-u(q) \ < \p-q\, V p,qe A ABC , 

and equality holds if and only if both p and q lie on a common segment £ a . This completes 
the proof. □ 
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As in [HI E2] , one can show by Lemma |5] that there is a unique measure preserving 
map T from (/, Aabc) to (g, Aabc) such that T (p) and p lie in a common £ a for all 
p G A^bc; an d satisfies the monotonicity condition 

(t (p) - r (g)) • (p - ?) > o Vp,qe£ a . 

With the help of and Lemma [U we prove that this T is indeed optimal. This fact is 
classical in optimal transport theory, but we show it in details for the sake of complete- 
ness. 

Lemma 7. To is an optimal mapping in the Monge mass transportation problem from 
(f,A AB c) to (g,A ABC ). 

Proof. Recall that the total cost functional is given by 

C(s)= [ f(z)\z-s(z)\dz, 



J Aabc 

where s G S, the set of measure preserving maps from (/, Aabc) to (g, Aabc)] an d the 
Kantorovich functional is defined as 



H^^)= g<p, 

J Aabc ''Aabc 

where (ip, ip) are function pairs in the set 

K, = {ip(x) +<p(y) < \x-y\ \/x,yeA AB c}- 
For all s G S and (if), tp) G K , we have 

(4.17) I{^ip) = j f(z)^(z)dz+ j f(z)<p(s(z))dz 

J A ABC J A ABC 

< [ f(z)\z-s(z)\dz 

J Aabc 

= C(s). 

That is 

sup / (ip, ip) < inf C (s) . 
k s 

Let u be the function constructed in the proof of Lemma [6], and let v = —u. Then we 
have (u,v) G /C. As T (p) and p lie on the same line segment, Lemma [6] implies 

u(z)-u (T (z)) = \z -T (z) \. 

So the inequality in (I4.17P becomes equality provided (if), <p) = (u, v) and s = T . There- 
fore 

C (T ) = I (u, v) < sup I (if), ip) < inf C (s) . 

K s 



REGULARITY IN MONGE'S MASS TRANSFER PROBLEM 



27 



Hence To is optimal and the segments £ a are transfer rays of Monge's problem. □ 

Let B' = (1, —4) be the reflection of the point B in the x-axis and let Q = Q* = Aabb 1 - 
Extend the functions /, g to fl such that they are symmetric with respect to the x-axis. 
From the proof of Lemma 5, one sees that /, g are smooth and satisfy the mass balance 
condition (jl.2p . The fact that r\ is quadratic close to shows that it can be reflected as 
a C 2 function, and Remark 1 shows that it is indeed smooth. It is also known (26] that 
Q and Q* are c-convex with respect to each other (for the cost function c e , < e < 1). 

Also extend T to Q so that it is symmetric with respect to the x-axis. By the 
uniqueness of monotone optimal mappings [T7], T is an optimal mapping of Monge's 
problem from (/, Q) to (g,£Y). 

We claim that To is not Lipschitz continuous at the point go = ( — 2, 0) . Let D a ^ be the 
strip in Aabc between the segments £ a and £ a +s, and let q a = (—2, a) be the intersection 
of £ a with the line {x = —2}, where 5, a > are constants. Let T (q a ) = (x a , y a ) . As T 
is measure preserving, we have (see the construction of the optimal mappings in [HI 132] ) 

lim - / f(x,y)dxdy = hm - / g(x,y)dxdy. 

6 JD atS n{x<-2} ° JD atS n{x<x„} 

That is 

/ (x + 2 + 3a)dx= [ (x + 2 + 3a)(l + ^-x + i](y/a(x + 2 + a))^dx. 

J-2-a J-2-a V 4 / 

Making the change t = 2 + a + x, we obtain 

ra rxa+2+a /l t — d \ 

/ (t + 2a)dt= / (t + 2a)[- + —— + rj (y/at) )dx 

Jo Jo \2 4 / 

Since both (t — a) and 77 (y/at) tend to when t, a — > (recall that rj(t) = 0(t 2 )), they 
are negligible in front of the constant |. This implies that, for small a, we should have 

(4.18) x a > -2 + (VE - 2)a. 

Indeed, either x a + 2 does not tend to 0, in which case ( 14. 18ft is satisfied, or it tends to 
0, in which case we can write, for small a, 



ra rXfj+2+a g 

/ (t + 2a)dt < / -(t + 2a)dx. 

Jo Jo 4 



Computing these integrals explicitly we get exactly the inequality (I4.18j) . 
On the other hand, by (]4.1j) . we have a = a 3 / 2 . Note that x(0) = —2. Hence 



x{a) - X M > I lim a -v 2 

cr^0+ 

Our claim follows. 



(4.19) lim " V " ; > - lim a~ 1 ' 2 = 00. 

0-^0+ a 4 <t->o+ 
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As go = ( — 2,0) is an interior point of Aabb>, we have thus constructed positive, 
smooth densities /, g, and c-convex domains Q = Q* = A^bb 1 , such that the associated 
optimal mapping T is not Lipschitz at interior points. 

As the triangle A^bb' is c-convex with respect to each other, the optimal mapping 
T e is smooth [26]. By Lemma 4, one has To = lim £ ^oT £ , and the above example shows 
that T £ is not locally, uniformly Lipschitz continuous as e — > 0. 

5. Applications and perspectives 

The regularity problem for the Monge cost is very natural in transport theory and 
very difficult. For the moment, even the implication f,g G C°° =r- T G C° in a convex 
domain is completely open. The transport To, among all the optimal transports for the 
cost \x — y\ (for which there is no uniqueness), is likely to be the most regular and the 
easiest to approximate. 

The present paper presented a strategy inspired by the previous results introduced 
in [26] to get Lipschitz bounds, i.e. L°° bounds on the Jacobian. Yet, it only allows 
for some partial bounds, and the counter-example of Section 4 shows that a Lipschitz 
result is not possible. However, in the same counter-example, the monotonic transport 
To is a continuous map, and the point where a non-Lipschitz behavior is observed shows 
anyway the behavior of a C°'5 map. Thus, it is still possible to hope for continuous, or 
even Holder, regularity results on T . 

We stress that these results could also be applied to the regularity of the transport 
density. The transport density is a notion which is specifically associated to the transport 
problem for the Monge cost (see p2]): it is a measure a which satisfies 

{div ■ (aDu) — f — g in f2 
\Du\ < 1 in Q, 

\Du\ = 1 a.e. on a > 0, 

together with the Kantorovich potential u. 

Several weak regularity results have been established, starting from the absolute con- 
tinuity of a if either / or g are absolutely continuous, till the LP estimates /, g G LP =>- 
a E LP (see pH Q3J O EH]). 

An explicit formula for a in terms of optimal transport plans or maps is available (we 
will not develop it here, see [I]) and most possible strategies for the regularity of the 
transport density need some continuity of the corresponding optimal transport. Yet, one 
of the advantages of working with a is that any optimal transport T produces the same 
density a. This allows for choosing the most regular one, for instance T , but requires 
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anyway some regularity on it. Here is where our analysis comes into play (without, 
unfortunately, providing any exploitable result). But there are other features of the 
transport density that one could take advantage of: from the fact that it only depends 
on the difference f — g, one can decide to add any common density to both measures. 
For instance, if / and g are smooth densities with compact support on M n , it is always 
possible to add common background measure on a same convex domain Q including 
both the supports. Q can be chosen as smooth as we want, and we can for instance take 
Q to be a ball. Also, one can add another common density to / and g so as to get g — 1. 
This last trick allows to avoid some of the tedious computations of Section 3, since in 
this case g(T) has not to be differentiated. 

In any case, even with these simplifications, the continuity result is not available for 
the moments. Possible perspectives of the current research involve the use of these 
partial estimates to prove continuity. 

Among the possible strategies 

• Use the bounds on DT e to get estimates on the directions of the transport rays 
for the limit problem, and use them to estimate how much the disintegrations of 
/ and g vary according to the rays. Using the fact that the monotonic optimal 
transport (in one dimension) continuously depends on the measures, one can 
hope for the continuity of T . 

• Use the fact that the bound on W gives an L°° bound on div(LDu) and, since 
L depends on \Du\, one faces a highly non-linear and highly degenerate elliptic 
PDE where the goal would be to get uniform continuity results on LDu. This 
recalls what has been recently done in very degenerate elliptic PDEs for traf- 
fic applications (see [29j K]); but seems (much) harder because LDu is not a 
uniformly continuous function of Du. 

• Write down some elliptic PDEs solved by some scalar quantities associated to 
T e , for instance by L, and use the bounds on the matrices A and w that have 
been proven here in order to apply De Giorgi-Moser arguments (or their wider 
generalizations, see [15] for a complete framework). Should it work, this would 
give Holder continuity. Unfortunately, our attempts have not given any useful 
PDE so far. 

All in all, up to the two-dimensional result of [19] (which requires disjoint and convex 
supports), the search for continuous optimal transports for the original cost of Monge is 
still widely open. 
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